Dark matter interpretation of recent electron and positron data 
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We analyze the recently released Fermi-LAT data on the sum of electrons and positrons. Com- 
pared to a conventional, pre-Fermi, background model, a surprising excess in the several hundred 
GeV range is found and here we analyze it in terms of dark matter models. We also compare with 
newly published results from PAMELA and HESS, and find models giving very good fits to these 
data sets as well. If this dark matter interpretation is correct, we also predict the possibility of a 
sharp break in the diffuse gamma ray spectrum coming from final state radiation. 



INTRODUCTION 

Recently, the question of dark matter in the universe 
has been put into focus, as data from a new generation of 
cosmic ray detectors have been tentatively interpreted in 
terms of indirect signatures of dark matter annihilation 
(or decay). In particular, although the data on antipro- 
tons from the PAMELA satellite [T] , and the mid-latitude 
data on diffuse gamma rays from Fermi show no hints 
of an "exotic", extra component, the positron fraction 
measured by PAMELA [5] and the sum of electrons and 
positrons by ATIC [4] definitely show an excess compared 
to standard background predictions. At any rate, these 
indications of an extra component has sparked a remark- 
able number of theoretical papers, both concerning dark 
matter explanations (for recent reviews, see [3 |6]) and 
other astrophysical sources such as pulsars (e.g., [7]). 

With the very recent release of Fermi-LAT [8| data on 
the sum of electrons and positrons ^ , the experimental 
precision has been increased in an impressive way. Al- 
though these high-quality data do not confirm the peak 
claimed by ATIC, they indeed also indicate an excess 
compared to the expected background as predicted in 
the so-called conventional, pre-Fermi GALPROP model 
[lOj . In this paper we show, however, that the dark mat- 
ter annihilation interpretation seems to fit even better 
with the new data. We compute confidence regions of 
the parameters needed in the dark matter sector to fit 
the data, and we show examples of the remarkably high 
quality fits obtained by some of the best-fit models. In- 
terestingly, the models that fit Fermi also reproduce the 
break indicated at high energy by the HESS collaboration 
[11 j . One should keep in mind though, that the pulsar 
interpretation seems still to be a possibility worthy of 
further studies [T5]. 

We point out two features of viable models, which had 
already been inferred in some other work by trying to 
fit previous data and checking other bounds (in partic- 
ular, from gamma-ray [T3J [TU [TS] and radio emission 
[m HZl HI] ) . Of course, the previous work tried to fit the 
"bump" claimed by the ATIC experiment, which is not 
seen in the new Fermi data. This means that smaller 
boost factors are actually required, and also the best 



fit mass may be increased. First, our solutions for the 
dark matter mass are in the 1-4 TeV range, which is 
indeed where the needed cross section enhancement by 
the Sommerfeld factor is operative [12] • Another mecha- 
nism giving a boost is the possible non-uniformity of the 
dark matter distribution within the nearest kpc of the 
solar system, or a combination of both, see, e.g., pOllST] . 
where possible boost factors much higher than needed in 
this work are obtained. 

Second, we find that the measured distribution seems 
to favour direct dark matter annihilation to with 
very little of e^e~ or t'^t~ , thus pointing to an unusual 
underlying particle physics model. However, if the anni- 
hilation is first into low-mass states which subsequently 
decay into leptons, both the absence of t+t~ (for kine- 
matical reasons) and the occurrence of /i^/i^ can be nat- 
urally explained, e.g. if the intermediate particle is scalar 
or pseudoscalar, by helicity arguments. 



DARK MATTER AND PROPAGATION MODELS 

We will consider two different types of dark matter 
models: i) leptonic models with direct annihilation to 
leptons, and as we will see, preferably to fi'^fi~, and ii) 
models annihilating to a light scalar (or pseudo-scalar 
or vector) decaying predominantly to leptons, where we 
focus on the Arkani-Hamed et al. (AH) type of models 
P2] and the Nomura and Thaler (N) type of models [15] . 
We use DarkSUSY [24] to calculate the positron fluxes for 
these models and for the case of the AH type models, we 
have modified DarkSUSY to allow for the scalar (vector) 
decaying directly to e~^e~. The propagation model we 
use is the one from [25| where we use the propagation 
parameters of the MED model of [26]. One should note 
though, that at the highest energies, which is of most 
interest in fitting the Fermi and HESS [TTl data, only 
the energy loss parameter tq plays a significant role and 
we choose tq — 10^^ s/GeV^ as a default value. For the 
dark matter halo model, to be conservative, we choose 
an isothermal sphere with a local dark matter density of 
Pa = 0.3 GeV/cm'^. However, our results are not very 
sensitive to the actual halo model, apart from the local 
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density po- For all classes of models, we have set the 
annihilation cross section to (crv) = 3 x 10~^^ cm^ s~^. 
Note that all the fluxes scale with Pq and at high ener- 
gies (all energies of interest here), they also scale with 
tq. With a boost factor from substructure and/or Som- 
merfeld enhancements. Bp, we can then define an overall 
enhancement factor, 

" (0.3 GeVcm-3) (1OI6 sGeV-O "-^^ 

which to a very good approximation will scale our fluxes 
up (or down). The enhancement factors we find below 
thus includes the Sommerfeld enhancement arising in the 
AH and N type of models. For both these latter models, 
we use the benchmark models AH1-AH4 and N1-N5 as 
defined in [TH] . 

We will also discuss other constraints, like constraints 
from gamma rays |13} I14[ 115] and synchrotron radiation 

nnnniiTH]. 

DATA AND BACKGROUNDS 

In our analysis, we will include the recent data from 
Fermi [3] and HESS [TT] for the summed spectrum of elec- 
trons and positrons, and Pamela data on the positron 
fraction [1]. We will not include ATIC |4] in our fits, 
as those data do not seem compatible with Fermi (al- 
though an excess is indeed also found by Fermi, the peak 
structure seen by ATIC is not). For HESS, we consider 
energies up to 2 TeV. Since the background model we 
use agrees with the HESS data (below 2 TeV) within 2cr, 
we use HESS data only to place an upper limit on our 
models, by requiring that our signal -I- background does 
not overproduce electrons at the HESS energies. 

For the background estimate, we use the conventional 
GALPROP [To] model (i.e. not specifically adjusted to 
fit the recent data; see also further discussion in the next 
Section), as given in 

As both data and background estimates have system- 
atic uncertainties, we have some freedom of adjusting 
their relative normalizations. It is a non-trivial task to 
do this in a fully consistent way, but looking at the back- 
ground estimate and the data it seems that the HESS 
data and the background estimate is slightly higher than 
the Fermi data. To take this normalization uncertainty 
into account, we have arbitrarily rescaled the HESS data 
and the background estimate by 0.85, which is within the 
expected systematic uncertainties of Fermi, HESS and 
the background estimates. 

RESULTS OF FITS 

We are now ready to perform our dark matter fits to 
the data. First we consider the leptonic models. It turns 



out that annihilation directly to e+e" typically gives a 
too pronounced peak at the WIMP mass Mdm- Tau 
leptons typically give a too soft spectrum and too many 
gamma rays from tt" decays from hadronic tau decays 
[131 [T7]. Hence, we will for these leptonic models focus 
on models annihilating only to It can be argued 

that this is unnatural (and not expected) for most dark 
matter candidates, but will come back to this question 
when we discuss the Nomura-Thaler and Arkani-Hamed 
type of models. 

In Fig. [TJl we show the required enhancement factor 
versus mass assuming 100% annihilation to /i^/i^. In[l^, 
the contours are shown for fits to Fermi and PAMELA 
data within 2a and the exclusion region which is 2a above 
the HESS data. As can be clearly seen, we typically get 
rather large regions with good fits to the data. The boost 
factors are of the order of 1000 and typically, slightly 
larger boost factors would be needed to get good fits 
also to the PAMELA data on the positron fraction. The 
quality of the fit is shown in Fig. [2^, where we show the 
spectrum for an example of a good fit model (the dot in 
Fig. [1^)). As can be seen, the excess in the Fermi data 
is very well-fit with this dark matter model and we get 
a reasonable (but a bit low) fit to the PAMELA data as 
well. [30! 

For the leptonic model, it is rather unnatural to as- 
sume that annihilation only goes to /i^/i^ . However, for 
the AH [22] and N [33] type of models, this follows nat- 
urally from the mass of the scalar being above 2m^, but 
below 2mr- In the N models, the light scalar (actually a 
pseudo-scalar) decays only to fJ,'^p~ , whereas in the AH 
models, the decay can go either purely to jjl^ jJC or to a 
mixture of e+e~ and fi^/i^ , depending on the nature of 
the scalar. Not surprisingly, the N models and the AH4 
model (which has a scalar that decays 100% to fJ.~^fJ.~) 
give the best fit to the electron and positron data and we 
therefore focus on these. In Figs.[T]3 and[2]3, we show the 
2a contours and an example spectrum for a good fit N3 
type of model where we have only kept the enhancement 
factor Ep and WIMP mass Mdm as free parameters (all 
the other N1-N4 models considered in [18] give the same 
positron fluxes and are thus equivalent in this respect). 

In Figs. [1]: and[2|:, we show the same 2a contours and 
example spectrum for the AH4 model. Even if the fits to 
the data for model N3 and AH4 are not as excellent as 
for the muon channel model, they are still good, and of 
course it is more natural to get predominantly muons in 
these models. Another advantage of these types of mod- 
els is that it is natural to get the required boost factors 
which enter the total enhancement factor (Eq. (jlj) from 
the Sommerfeld enhancement. 

Of course, one has to investigate other constraints on 
these kinds of models, especially given the rather large 
boost factors needed to fit the electron and positron data. 
For the leptonic model, we can use the constraints in [IB], 
where the gamma ray and radio constraints, primarily 
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FIG. 1: The 2a contours in the enhancement factor - mass plane for a) annihilation to fi'^ fi~ , b) the Nomura- Thaler model N3 
and c) the Arkani-Hamed et al. model AH4. The contours are shown for PAMELA and Fermi, whereas the HESS data is only 
used as an upper limit. The black dot is the example model shown in Fig[2] 
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FIG. 2: Spectra for examples of good fit models in[l The signal and background are shown for electrons (e"*" + e~) together 
with Fermi [5] and HESS data [111 I27| . The HESS data and the background model has been rescaled with a factor 0.85. In 
the inset, the positron fraction as measured with PAMELA is shown together with the predicted signal for the same model. 



towards the galactic centre and dwarf spheroidals were 
investigated. For Einasto or NFW profiles, the best fit 
models are excluded due to gamma rays from the galactic 
centre. However, for less steep profiles, like an isothermal 
sphere, our best fit models are not excluded by these 
data. 

For the N and AH models, constraints from gamma 
rays and radio (including final state radiation photons) 
were investigated in [18]. The same conclusion holds for 
these models, if the halo profile is an Einasto or NFW 
profile (or steeper), the models are already excluded. 
However, for shallower halo profiles, like an isothermal 
sphere, the models are still viable. One should note that 
the electron and positron fluxes discussed in this paper 
are not very dependent on the choice of halo profile, so 
the best-fit models derived here, would be more or less 
the same for an NFW profile instead of the isothermal 
profile we used in our analysis. 

Given the large amounts of high-energy electrons and 
positrons injected into the galaxy with these models, it 
is also fair to wonder about secondary radiation from 
inverse Compton scattering on the interstellar radiation 
field mimilTlllH]. In [g it is concluded that models 
annihilating to n^fJ,~ are at tension with EGRET data 



and that Fermi will be able to probe these models. Given 
the new Fermi data, lower boost factors are needed than 
those assumed in [2] , so the tension with EGRET data is 
less severe. However, Fermi should still be able to probe 
these models. For the N3 and AH4 model, we get very 
similar constraints [T7] and these are also viable with a 
shallow halo profile. 

One should also note that we have chosen to work 
with a rather standard halo and diffusion model, but it is 
rather straightforward to rescale our results via the en- 
hancement factor introduced in Eq. ([l]). Note that the 
dependence on po and tq in Eq. ([T|) is a very good ap- 
proximation for high energies. For lower energies (i.e. the 
PAMELA range), it is more involved as the positrons at 
these energies have propagated rather far. Keeping the 
signal fixed at higher energies, it is possible to move the 
signal from dark matter up at lower energies by having 
a larger significant diffusion region (by having a larger 
diffusion zone half height and a larger diffusion coeffi- 
cient). Increasing tq will also increase the ffuxes at low 
energies slightly more than the linear relation in Eq. ([T]) 
as positrons then sample a larger (and partly denser) re- 
gion in the galaxy. These effects are more pronounced 
for steeper halo profiles, like a Navarro-Frenk- White [59] 
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profile, than the isothermal profile considered here. In 
either case, these effects are rather small, and Eq. ([T]) re- 
mains a very good approximation also for the PAMELA 
energies, even for steeper halo profiles. 

CONCLUSIONS AND DISCUSSION 

We see that it is certainly possible to fit the excess 
seen in the Fermi data with a dark matter model. We 
have here used a 'standard' diffusion model and a 'stan- 
dard' dark matter halo, but one should keep in mind that 
we need to invoke boost factors (either from substruc- 
ture and/or a Sommerfeld enhanced annihilation cross 
section) of the order of 10'^ for this to work. For the 
dark matter models annihilating directly to leptons, we 
also prefer models annihilating predominantly to 
as the e^e'" channel gives too many hard electrons and 
the r+r" channel gives too soft electrons (and possibly 
too many gamma rays). Hence, for these models, we 
would need to invoke a rather unnatural selection of the 
/Lt^/x" channel. However, for the Sommerfeld enhanced 
models with light scalars, we naturally get annihilation to 
only (or at least predominantly) fJ.^ fJ.^ through annihila- 
tion to light scalars which then (for kinematical reasons) 
decay mostly to muons. These models typically require 
a slightly larger mass and enhancement factor, but on 
the other hand in these models the enhancement factor 
from Sommerfeld enhancement is expected, a we get the 
annihilation to muons 'for free'. 

One should also note that we have not varied the halo 
or propagation model parameters, which could help in 
getting an even better fit also to PAMELA data, as dis- 
cussed above. 

We have here used the conventional, pre- Fermi, con- 
ventional GALPROP [To] model for our background esti- 
mate. This model should most certainly be revised given 
the new data, but we leave this for future work. 

A very nice property of DM models, as was recently 
pointed out |21j and as verified in this work, is the relative 
model independence of the shape of the DM spectrum at 
low energies. The case for a dark matter interpretation 
of the new surprising data is in our view getting stronger, 
although we have to await even more data, plausibly us- 
ing gamma-rays, to settle the issue. 

As a final note, for the very best fit model that an- 
nihilates directly into muons (or, but to a lower extent, 
models with spin-0 bosons decaying to muons) , there may 
in fact exist a unmistakeable signature that would make 
it uniquely different from the pulsar interpretation (or 
other backgrounds), namely the direct emission of final 
state radiation (FSR) in the annihilation process. Al- 
though the authors of ^4j unfortunately did not add that 
contribution to the total flux in their Fig. 1, they show 
its magnitude as a dashed red line. It is obvious that a 
very striking feature should appear - a slow rise of an 



excess in the E'^dN^/dE curve up to 100-300 GeV where 
the rise should suddenly be much steeper as the FSR 
dN^/dE (X 1/E sets in, i.e. a linear rise when multiplied 
by E'^ (see Fig. 1 in [TB]). However, the existence and 
exact location of this feature will depend on the rather 
poorly known starlight energy density near the galactic 
centre, and also the detailed spatial distribution of the 
DM halo. It will be very interesting to see if Fermi-LAT 
can find (or rule out) such a behavior when they analyze 
the diffuse gamma radiation from the region nearer the 
galactic centre in more detail. At mid-latitudes and at 
lower energies such an excess is neither present [2] nor 
predicted to exist in the models discussed here. 
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